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Abstract 



The exact Quantum Kinetic Equations for a two-flavour active-sterile neu- 
trino system are used to provide a systematic derivation of approximate evolu- 
tion equations for the relic neutrino asymmetry. An extension of the adiabatic 
approximation for matter-affected neutrino oscillations is developed which in- 
corporates decoherence due to collisions. Exact and approximate expressions 
for the decoherence and repopulation functions are discussed. A first pass is 
made over the exact treatment of multi-flavour partially incoherent oscilla- 
tions. 



1 



I. INTRODUCTION 



Taken together, the solar neutrino |IJ, atmospheric neutrino Pj and LSND J| results 
require at least four light neutrino flavours if neutrino oscillations are to be their explana- 
tion. Consistency with the measured invisible width of the Z boson then demands that 
any additional light flavours be sterile with respect to ordinary weak interactions. With 
these important experimental results as partial motivation, the cosmological consequences 
of sterile neutrinos have recently been re-examined f4|-|TT||. It has been shown H that os- 
cillations between active and sterile neutrinos can create significant asymmetries between 
the number densities of relic neutrinos and antineutrinos.Q These asymmetries then have 
extremely important consequences for Big Bang Nucleosynthesis (BBN), through their sub- 
sequent suppression of active to sterile neutrino oscillations ]7[H|, and through the effect on 
BBN nuclear reaction rates of a large v e asymmetry [P,|10|l. 

Various approximations have been explored in the study of neutrino oscillations in the 
early universe, with neutrino asymmetry or lepton number evolution being an important as- 
pect of these investigations |T^-[14| . The pioneering analyses neglected the thermal spread in 
relic neutrino momenta: the evolution of the neutrino ensemble was assumed to closely track 
the evolution of neutrinos having the mean momentum (p) ~ 3.15T. Other approximations 
often utilised were the rate equation approximation also neglecting the thermal momentum 
distribution, or a Pauli-Boltzmann approach when the momentum spread was not neglected. 
In the case of active-sterile oscillations, the number densities of sterile neutrinos were often 
neglected if they were small. Pauli blocking, except for the BBN reactions, has usually been 
neglected. Various other more subtle approximations have also been made. For instance, in 
the presence of a nonzero neutrino asymmetry, the collision rate for neutrinos differs from 
that of antineutrinos. This difference has always been neglected. The repopulation of active 
neutrino distributions by weak interactions is another process often treated approximately. 

The main purpose of this paper is to clarify the nature of these various approximations, 
and to discuss their regions of applicability more fully and systematically than has been done 
in the past. Our approach will take the known exact Quantum Kinetic Equations (QKEs) 
ni.[TH] for a two-flavour active-sterile system and then develop systematic approximations 



to them, focussing on the evolution of lepton number. Approximate treatments of lepton 
number evolution have in the past been vital in obtaining analytical insight. In fact, the 
discovery that partially incoherent active-sterile oscillations can create lepton number expo- 
nentially quickly was made through a treatment that neglected the momentum spread and 
treated the evolution in the centre of the MSW resonance very approximately |J . The effect 
of removing these approximations was subsequently investigated, and a more exact treat- 
ment of lepton number creation obtained |7J]. Nevertheless, the simple and approximate 
starting point proved essential in the discovery of this interesting effect, whose existence 
had hitherto been overlooked because of the complicated nature of more exact treatments. 
Another important motivation for studying approximation schemes is the computationally 



1 Synonymns for "relic neutrino asymmetry" are: (i) neutrino chemical potential (when thermal 
equilibrium holds), (ii) neutrino degeneracy, and (iii) nonzero lepton number for the universe. 
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intensive nature of numerical analyses based on the exact QKEs. The neutrino momentum 
spread is computationally demanding, especially during the BBN epoch when neutrinos are 
out of thermal equilibrium. 

Another goal of this paper is to address multi-flavour effects in the early universe from 
first principles for the first time. Previous analyses have either considered a two-flavour 
system in isolation, or have approximated the multi-flavour system via pairwise two-flavour 
subsystems |7],^,|n]]. (Note that some genuine three-flavour effects were considered for fully 



coherent oscillations in Ref. ||10|| .) We will clarify how the "pairwise two-flavour approx- 
imation" arises from a more exact treatment. The three-flavour problem is particularly 
important, because one of the principal applications of lepton number creation lies in the 
suppression of the large angle — v s oscillations that may solve the atmospheric neutrino 
anomaly ||16|| . The lepton number responsible for this suppression is quite possibly or plau- 
sibly created by small angle v T — v s oscillations |7|||. If so, then a full three-flavour analysis 
of the Up, v T and v s system is warranted. 

This paper is structured as follows: Section discusses exact and approximate treat- 
ments of the evolution of lepton number in the context of a simple two-flavour active-sterile 
system. Section |T| briefly discusses the important case where the strictly sterile neutrino is 
replaced by a mirror neutrino. Section [TV] is devoted to a study of multi-flavour effects, and 
Section [V] is a conclusion. 



II. THE TWO-FLAVOUR ACTIVE-STERILE NEUTRINO SYSTEM: FROM 
FIRST PRINCIPLES TO USEFUL APPROXIMATIONS 

The early universe contains, in part, an ensemble of neutrinos.^ Their evolution in time 
is affected by three phenomena. The first is simply the expansion of the universe. The 
second is decohering collisions with the background medium. The third is, of course, the 
coherent process of oscillations governed by a matter-dependent effective Hamiltonian. The 
effect of collisions and oscillations is quantified using the Quantum Kinetic Equations. The 
QKEs generalise the Pauli-Boltzmann Equations to include quantum coherence between 
the particle species, in our case the active and sterile neutrinos (and active and sterile 
antineutrinos). 

We will focus on the two-flavour system comprising one active and one sterile neutrino 
species in this section. We do so because two-flavour active-sterile oscillations lead to the 
rapid creation of lepton number when the oscillation parameters are in the correct region. 
The multi-flavour case is deferred to Section |V]. Our task is to track the evolution of the 
neutrino and antineutrino ensemble, focusing in particular on the evolution of lepton number. 
We will take as our initial, high temperature, condition that the number densities of sterile 
neutrinos and antineutrinos are zero. This is expected simply because, by definition, sterile 



2 When there is no chance of confusion, or when the distinction is unimportant, we will use the 
term "neutrino" to mean neutrino and/or antineutrino, active and/or sterile. When the distinction 
is important, we will specify which one we mean. 
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neutrinos do not feel electroweak (or strong!) interactions and so will decouple very earlyfj 
We will focus on the epoch between T = m M ~ 100 MeV and BBN (T ~ 0.1 - 1 MeV), 
so that the background plasma consists of essentially only photons, electrons, positrons, 
neutrinos, antineutrinos and whatever amount of sterile neutrinos and antineutrinos get 
generated through oscillations. 



A. The exact Quantum Kinetic Equations 



We first write down and discuss the exact QKEs, partly by way of review and to de- 
fine notation, and partly to discuss approximations to the decoherence and repopulation 
functions (to be defined shortly). 

Consider a two-flavour active-sterile system composed of u a and u s , where a denotes 
either e, p or r. The 1-body reduced density matrices describing v a — v s and v a — V s 
oscillations are given, respectively, by 

P (p) = ip (p)[i + P(p)-ff], (i) 
p(p) = ip (p)[i + F(p)-*], (2) 

where 



P(p) = P x (p)x + P y (p)y + P z {p)z, (3) 

plus an analogous expression for antineutrinos. The functions Pq and P are just the coeffi- 
cients of a convenient expansion of p in terms of Pauli matrices a and the identity matrix (P 
is often called the "polarisation"). The quantity p is the neutrino or antineutrino momen- 
tum. It will be understood that Pq and P depend on time (or, equivalently, temperature via 
t ~ mpi anc k/HT 2 ) as well as momentum. The diagonal entries of p (p) are relative number 
density distributions in momentum space of u a (z7 Q ) and v s (z7 s ): 

N a (p) = ^P (p)[l + P z (p)]N^(p,0), (4) 
N s (p) = l -P,{p)\l - P z (p)}N eq ( P ,0), (5) 

where N eq (p, 0) is a reference number density distribution that we have set equal to the 
equilibrium Fermi-Dirac distribution, 



3 Negligible initial sterile neutrino number densities is a consistent assumption because the neutrino 
system is frozen at high temperatures — see later. It is at lower temperatures that sterile neutrino 
production due to partially incoherent active-sterile oscillations threaten overproduction. 
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with zero chemical potential. This choice amounts to normalising p such that Tr(p) = 2 
if thermal equilibrium holds and the chemical potential is zero. The antineutrino relative 
number densities are given similarly. Since the physical interpretation of p is related to 
number density distribution ratios, the expansion of the universe does not directly contribute 
to the time evolution of p. Note that N eq (p, 0) depends on time only through the expansion 
of the universe, and not on coherent or incoherent matter effects. 
The evolution of Po(p) and P(p) is governed by the QKEs 



^P(p) = V(p) x P(p) + [1 - P x (p)\ 



|lnP (p) 



£>(p) + ^lnP (p) 



[P x (p)x + P y (p)y) 



and 



^Po(p) = Rip). 



(?) 



(8) 



These equations are obtained by evolving the full density matrix for all particles in the 
plasma forward in time by means of the S matrix, and then tracing over all degrees of 
freedom other than u a and v s . The essential difference between the Quantum Kinetic and 
Pauli-Boltzmann approaches is that the former time evolves amplitudes while the latter 
evolves probabilities. See, for example, Ref. []14| for a detailed derivation of the QKEs and 
further discussion. The antineutrino QKEs take an identical form, but with the substitutions 



P, V^V, D^D, R^R. 



(9) 



We now discuss the form and significance of each of the terms above.0 

The function V(p) describes the quantally coherent part of the matter- affected evolution 
of the v a — v s subsystem. It is given by 



= f3(p)x + \(p)z, 



(10) 



with 



Pip) 

Hp) 



Am 2 
2p 

Am 2 



2p 



sin 29 , 
cos 2^ + V a . 



(11) 
(12) 



The quantities Am 2 and 8q are the mass-squared difference and vacuum mixing angle for 
v a — v B oscillations, respectively. We define the mass eigenstate neutrinos v a h by 



4 These equations are exact for strictly sterile species. All contributions due to possible right- 
handed weak interactions and the like have been neglected. 
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u a = cos 9 v a + sin 9 u b , v s = — sin 9 i> a + cos 9 u b , (13) 

with #o denned so that cos2#o > and Am 2 = m 2 — m 2 . The function V a is the effective 
matter potential |17|]. To leading order, including the lowest order finite temperature term, 
it is given by [JTEf 

K, = -^-[-a(p)+6(p)], (14) 



with the dimensionless functions a(p) and 6(p) given by 

a(p) = 



4C(3)V2G F r 3 p r(Q) 

tt' Am* L ' (15) 
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4C(3)y/2A a GpT A p 

Hp) = 2 A 2 2 » ( 16 ) 

where T is temperature, is the Fermi constant, my/ is the VT-boson mass, £ is the 
Riemann zeta function, A e ~ 17, ~ 4.9 and 

= L a + L e + L^ + L T + V . (17) 

The quantities on the righthand side of this equation are the various asymmetries, given by 

n t — ri-j / x 



where 



n 



f 



N f (p)dp (19) 



are number densities (note that we will abbreviate L Ua by L a in this Section). The a{p) 
term in Eq. (|T4"D is the leading order density-dependent (Wolfenstein) contribution to the 
effective potential, while the b(p) term is the leading order finite temperature term. (Note 
that the Wolfenstein term also depends on temperature, but only indirectly through the 
cosmological red-shifting of fermion number density distributions.) Observe that V depends 
on p through the dependence of a(p) on L a . As we discuss later, this makes the evolution 
of L a non-linear. It is important to notice that the dependence of V on L a is an 0(Gp) 
effect (due to coherent forward scattering), rather than an 0{G 2 F ) effect. The quantity rj is 
a small term related to the cosmological baryon-antibaryon asymmetry. For antineutrinos, 
the corresponding function V is obtained from V by setting 

A(p) = -^[cos2^ - b(p) - a(p)], (20) 

and replacing A by A. In other words, one simply replaces by —L^ a \ 

The Mikheyev-Smirnov- Wolfenstein (MSW) [0,11^] resonance conditions are 

cos 2# + a(p) - 6(p) = 0, cos 26 - a(p) - b{p) = 0, (21) 
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for neutrinos and antineutrinos, respectively. Note that prior to the creation of lepton 
number, both resonance conditions are simply b(p) = cos2# , which can only be satisfied if 



Am 2 < 0. (22) 

This condition plays an important role (see later). It is also important to realise that the 
resonance conditions of Eq.(^TJ) at a given temperature T are met only for neutrinos and 
antineutrinos of a certain momentum p = p res . This is why the spread of neutrino momenta 
plays a significant role in the analysis M. 

The decoherence or damping function D{p) quantifies the loss of quantal coherence due 
to u a collisions with the background medium. Its exact expression is derived in a general 
form in Ref. ||14| , and is given in the Appendix by Eq. ( |A1| ). When thermal equilibrium 
holds, the general expression for D reduces to a useful compact form, given by 

D{p) = \T a {p) (23) 

where T a {p) is the total collision rate for i/ a 's with momentum p. By examining the mo- 
mentum dependence of the specific weak collision processes operating between the neutrino 
decoupling temperature and 100 MeV, we see that 

r a ( P ) = (r a )JL, (24) 

where 

< P> °=18&) T " 3 ' 15T <25) 

is the average momentum for a relativistic Fermi-Dirac distribution with zero chemical po- 
tential (indicated by the subscript "0"), and (r a ) is the thermally averaged total collision 
rate for v a . This last quantity can be expanded as a power series in L a , 

(r«) = G 2 F T 5 y a (l - z a L a ) + 0(L 2 a ), (26) 

with y e ~ 4, y^ T ~ 2.9, z e ~ 0.1 and z^ T ~ 0.04. The y a term, 

(r a ) = y a G 2 F T 5 , (27) 

is the thermally averaged collision rate at zero chemical potential, while the z a term is the 
first order correction due to a finite neutrino asymmetry. For antineutrinos, the correspond- 
ing expression is 

Dip) = ±T a ip) = l -G 2 F T"y a ^il + z a L a ) + 0{L 2 a ), . (28) 

In previous studies Eq. (p4|) was adopted as an approximate result for T a (p). Actually, 
it turns out to be exact, as we show in the Appendix. 
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It is important to observe that although D is in principle a dynamical quantity through 



its dependence on neutrino number densities and hence on p [see Eq.(Al)], it becomes ap- 
proximately kinematic provided two conditions hold: (i) thermal equilibrium and (ii) lepton 
number being small. If these are so then D(p) — y a G 2 F T 5 p/2(p) . The fact that D is most of 
the time just a given external function of temperature and momentum rather than a dynam- 
ically evolving quantity simplifies the numerical solution to the equations. Observe that the 
absence of thermal equilibrium is correlated with sufficiently weak interaction rates, so D is 
necessarily small in that case relative to V (for relevant choices of oscillation parameters). 
If, on the other hand, thermal equilibrium holds, but the neutrino chemical potentials are 
nonzero and evolving in time, then the dependence of D on p will only be important when 
the asymmetries are very large. For interesting choices of oscillation parameters, this only 
occurs at temperatures approaching the BBN epoch when collisions are again becoming 
insignificant. It turns out that in practice one may neglect all but the leading neutrino 
asymmetry independent term y a G F T 5 p/2(p) in the applications we have in mind, a point 
we will discuss more fully later on. 

Pauli blocking is neglected in the calculation of D, and also the repopulation func- 
tion R discussed below. To incorporate this effect, appropriate factors of (1 — /), where 
f(p) = 2n 2 N(p)/p 2 , need to be inserted into the integrals. To estimate the size of the error 
introduced by neglecting Pauli blocking, we can evaluate / at the peak of the equilibrium 
momentum distribution, p ~ 2.2T. Since /(2.2T) ~ 0.1 we estimate the error to be at most 
10%. 

We now discuss the repopulation or refilling function R(p) which controls the evolution 
of Po(p). Since 

the evolution of Pq is governed by processes that deplete or enhance the abundance of 
j/qj's with momentum p. The rate of change of Pq(p) clearly receives no contribution from 
quantally coherent v a — v s oscillations, because the two flavours simply swap. It therefore 
equals the rate at which u a 's of momentum p are generated by scattering processes, minus 
the rate at which f a 's of momentum p are scattered out of that momentum value. Its general 
form is 

n/n 27r r d 3 k' d 3 p' d 3 p . . 

E{t' ! KW.iWI".(*').j(p')l[/».('=')/ J (!'') - 

3 

+ v 2 [v a (k),v a (p)\j(k% mUWjip') - UMfMh (so) 



where 



1 + exp( 



T 



and V(i, j\i' , j') denotes an interaction matrix element (defined in the Appendix) for the 
process % + j —^i' + f. 
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As for the decoherence function D, the general expression for R in Eq. (p0|) simplifies 
when thermal equilibrium holds. For temperatures between neutrino decoupling and 100 
MeV, it is given by 

R{P) = Ta{p) " \ P ^ )[l + P * {P)] } ' (32) 

where N eq (p, //) is the Fermi-Dirac momentum distribution function with chemical potential 
/i, and /Iq, is the chemical potential for v a . This expression follows from Eq.(^) when 
all number densities except for v a are thermal, and v a is approximately thermal (see the 
Appendix for the derivation). The chemical potential /i a is obtained through the lepton 
asymmetry using 



T 3 



6n 7 V T 



(33) 



The physical interpretation of this expression for R is that all weak interaction processes 
involving v a are tending to send its number density towards the equilibrium configuration. 



Equation (|32| ) was adopted in previous works on a heuristic basis ||10|| . In fact, it is an 
exact result under the conditions stated above (see the Appendix). Observe that above 
the neutrino decoupling temperature, v a is approximately in thermal equilibrium and so 
R is very small (this issue is discussed more carefully in the next Subsection). Below the 
neutrino decoupling temperature, the form given in Eq. (^) should in principle be used. For 
antineutrinos, R is obtained from R by replacing [i a by —\i a (when thermal equilibrium 
holds). 

Clearly, the dependence of the neutrino momentum distributions on the repopulation 
function will be most sensitive at low temperatures and when the asymmetry has evolved to 
large values. At high temperatures, the weak interaction rates are fast enough to thermalise 
the distributions rapidly, spreading out an asymmetry across the distribution (as described 
by the chemical potentials). At lower temperatures, this does not occur as efficiently, so 
the correct form for the repopulation function becomes more complicated [see Eq.(^)] if we 
want to track the actual form of the neutrino momentum distributions. We will discuss the 



low temperature (T ~ a few MeV) regime later on in Subsection HE 



B. Evolution of neutrino asymmetry 

We now focus on extracting the evolution of lepton number or neutrino asymmetry L a 
from the Quantum Kinetic Equations. Recall that this is a crucial quantity because the 
Wolfenstein a(p) term in the effective matter potential [see Eq.([14])] is proportional to a 
linear combination of fermion, and in particular neutrino, asymmetries. If this term is 
appropriately large, then active-sterile oscillations will be suppressed at low temperatures. 
Also recall that L e affects primordial light element abundances through weak interaction 
BBN reaction rates. 

In principle, the QKEs (|7|) and (|j) (plus the antineutrino equations) are solved to in- 
directly yield the evolution of the asymmetry defined by Eq. (|T%|) . It is, however, very 
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important to obtain simpler approximate evolution equations for L a . There are two reasons 
for this. First, it is difficult to understand the qualitative behaviour of the evolution of the 
asymmetry from the in-principle QKE procedure per se. Analytic insight requires that the 
evolution of L a be directly considered and approximations introduced. Second, the direct 
numerical solution of the full QKEs is computationally demanding because of the neutrino 
momentum spread and the MSW resonance phenomenon. 

We begin our study of approximate evolution equations by observing that Eqs.(|5]), (|1 
and (|19D combine to produce 



~dt 



dPo ( 
dt 1 

dt 



1 + P Z 



P 



dP z dP „ , - dP z 

ir-ir< 1 + p '> +/> °— 



J I3(P„P, - P„P„) 



2-^ + P(P P y - P Q P y ) 



dt 

N e \p, 0)dp 



N c %p, 0)dp 



n- 



ri- 



(34) 



as the time evolution equation for the a-like neutrino asymmetry L a . We have used the fact 
that N eq (p, 0)dp/n 7 is time independent, and also conservation of lepton number through 



N eq (p, 0)dp n a + n s — n a — n s 



constant. 



ri- 



ii-, 



(35) 



So far, no approximations have been made (except that possible tiny flavour changing neu- 
tral current effects for massive neutrinos have been neglected.^) Using Eq.(p9[) and the 
corresponding expression for antineutrinos, Eq.(|34"D implies that 



dL a 



2n„ 



(3 (N a + N s )P y - (N a + N s )P y dp. 



(36) 



To proceed, we need approximate solutions for P y and P y . These are obtained by solving 
the QKEs for P and the corresponding antineutrino equations. In expanded form, the QKEs 
(0) are 



dP x 
dt 

d_Py 

dt 
dt 



-XP y -DP x -^^. 
v x P dt 



XP X - (3P Z - DP V 



PydPo 
Pn dt ' 



0Py + 



P z dP 



Po dt 



(37) 



which also require Eq.(||). 

The qualitative character of the evolution goes through distinct phases. We first discuss 
the high temperature initial conditions. As stated previously, it is expected that the initial 



5 These effects will always be small, with their exact size being model dependent. 
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abundance of sterile species is very close to zero. This implies that P z ~ 1 at high T, so that 
N s ~ 0. At temperatures above the neutrino decoupling temperature, the neutrinos are by 
definition in thermal equilibrium, so the repopulation function R equals zero according to 
Eq . (p2|) . This means that Pq is a constant which depends on the initial neutrino chemical 
potential and hence on the initial neutrino asymmetry. We will also take the initial neutrino 
asymmetry to be negligible (it could be of the order of the baryon asymmetry, for instance^). 
If this is so, then Pq ~ 1. 

The first phase occurs at high temperatures. With the asymmetry set to its initial value 
of zero, the function A is dominated by the finite temperature term of the effective potential 

Am 2 , , , , , 

V a ~ — b (p . 38 
2p 

Replacing p by its thermal average, we see that the magnitude of this term increases with 
temperature as T 5 . The function (3 effectively goes to zero at sufficiently high temperature 
because it scales as T -1 . The decoherence function D, like A, scales as T 5 . (Note that, 
numerically, |V^/D| — 60 or 24 for a = e or fi/r, respectively.) Setting P to be constant, 
and neglecting f3, the QKEs ( pTf ) become 

BP BP BP 

^L~-\P y -DP x , ?±l~\P x -DP y , ^~0. (39) 
Bt y Bt * i "at ( J 

Therefore P z {t) is frozen at its initial value P z (t) — 1, and no sterile neutrinos are produced. 
Writing 

P x + iP y = \P\e l t (40) 
the approximate QKEs (|3~9D become 



B\P\ 



-D\P\, ^~A, (41) 



Bt 1 " Bt 

which are trivially solved to yield 

|P(t)|e^ (t) ~ |P(0)|e^°)e-/o ^^e^o x ^ dt ' . (42) 

Since, at high temperatures, D is large compared to the expansion rate of the universe we 
see that \P\ and hence also P x and P y get exponentially damped to zero very rapidly from 
any initial values. (Observe that they rapidly oscillate about zero as they do so due to 
the large A.) The system is therefore completely incoherent ("in a mixed state") at high 
temperatures, and frozen (Quantum Zeno or Turing Effect). 



6 A natural possibility is that the primordial neutrino asymmetry is produced by a similar mecha- 
nism to the baryon asymmetry. It is of course also possible that some relatively high temperature 
mechanism operates that produces a primordial neutrino asymmetry that is much larger than the 
baryon asymmetry. However, we will not consider this scenario in this paper. 
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The third phase occurs at very low temperatures where the collisional terms D and b 
are negligible. The neutrino subsystem then evolves as a coherently oscillating subsystem 
which is also coupled through R to the background. The oscillations are matter-affected if a 
sufficiently large neutrino asymmetry has been created by this time. This is the regime 
relevant for the BBN epoch, for interesting choices of oscillation parameters. We will discuss 
this regime in a later Subsection. 

Intermediate between these two limiting regimes lies the second phase. During this phase, 
the system emerges from its frozen initial state. The evolution is driven by an interplay 
between collisions, which are still rapid compared to the expansion rate, and oscillations, 
which begin to affect the neutrino number densities through a non-negligible (3 (recall that 
(3 is proportional to sin 2 26$ where 9q is the vacuum mixing angle). This is the regime of 
the "static approximation" introduced in Ref. J7|. It is during this epoch that neutrino 
asymmetries will evolve to values orders of magnitude higher than the baryon asymmetry 
provided the oscillation parameters are in the appropriate region. We will now provide new 
analytical insight into this phenomenon by deriving in a new way the very useful approximate 
form for the asymmetry evolution equation considered in Ref. 0. We will see that the 
static approximation is actually a combination of an adiabatic approximation for partially 
incoherent oscillations, and a small (3 expansion. This point was not crisply realised hitherto. 

We now take the first crucial step, that 



OPo 
dt 







continues to be a good approximation. With this approximation, the QKEs (|3 
the homogeneous equations 



(43) 
simplify to 



d_ 

dt 





' —D 


-A 




fP x 


(a)-( 




—D 


?) 






v 


(3 




u 



(44) 



or, in a self-evident matrix notation, 



dP 

~dt 



/CP. 



(45) 



The continued justification for Eq.fl4"3]) is as follows: Above the neutrino decoupling tem- 
perature, most of the v a ensemble is in thermal equilibrium. From Eq.(|32|) we see that R 
is therefore zero or very small. One might be tempted to conclude that in fact R is ex- 
actly zero above the neutrino decoupling temperature. This, however, cannot strictly be 
the case. Suppose u a 's of some given resonance momentum p res oscillate strongly to u s 's. 
Instantaneously, N a (p res ) goes to zero or close to it, but Po(Prcs) does not change due to 
these oscillations because of the generation of a nonzero N s (p ies ). However, the absent mo- 
mentum mode of z/ Q 's is quickly repopulated from the background medium, leading clearly 
to an overall nonzero value for ^ evaluated at the momentum p ies . We would therefore 
expect that R is very small above the neutrino decoupling temperature, except when the 



oscillations are very strong. This means that the approximation in Eq.(f4"3"|) will be a good 
one except possibly in the centre of an MSW resonance when oscillations might be rapid. 
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The word "might" is used because the oscillations will be non-adiabatic in some regions of 
parameter space. So, we proceed with the understanding that our approximation scheme 
might not be strictly valid in the centre of a resonance. We will have cause to revisit the 
centre of the resonance later on. 

To solve Eq.(H3), we first introduce the "instantaneous diagonal basis" through 



(Qi\ 
Q2 



Q = WP, 



(46) 



where U is a time-dependent matrix that diagonalises /C, 

JCd = diag(/ci, k 2 , fc 3 ) = WCW 1 , 
with ^1,2,3 being eigenvalues. In the instantaneous diagonal basis, Eq.flH]) becomes 



dQ 

dt 



JC d Q 



(47) 



(48) 



This equation resembles, but is not the same as, the MSW evolution equation written in the 
instantaneous mass eigenstate basis with a varying background density. In that case, the 
adiabatic approximation sees the time derivative of the instantaneous mixing matrix set to 
zero. 

Our next important approximation after Eq.(H3T) is to analogously set 



~dt 



— ~ 0, 



so that the evolution equation becomes 



dQ 

dt 



JC d Q. 



(49) 



(50) 



This is clearly within the family of adiabatic-like approximations. (We defer the discussion 
of the region of applicability of this approximation to a later subsection.) The difference 
with the usual adiabatic approximation for MSW evolution is the existence of quantum de- 
coherence through a nonzero D. We will now show that Eq.(|50|) reproduces the approximate 
evolution equation derived in Ref. @ using what was termed the static approximation. The 
static approximation of Ref. /fj/ is therefore revealed to be just the adiabatic approximation 
for a partially incoherent system of neutrinos undergoing matter affected oscillations (in the 
small f3 limit — see below). (Note that this type of idea was briefly discussed in Ref. [TBI]-) 
Equation (150) can be formally solved to yield 



Qi(t) 

Q2(t) 

Qs(t) 



( e f Q kv{t')dt' 



V 














J*k 3 (t')dt> 



Qi(o) 
Q 2 (o) 
Q 3 (o) 



(51) 



which in turn implies that 
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Py(t) 


= U~\t) 







( Xk^df 



\ 










e fi 
e Jo k ^') dt ' 





"Px(0)" 


W(0) 









(52) 



The mixing matrix W 1 is obtained by placing the normalised eigenvectors Ki in the columns, 
where 



1 

Mi 



1 



"0 



A 

D+kj 
\ki 



(53) 



and 



Alfcl 



A 2 |fci| 2 + (/3 2 + |A; i | 2 )|D + hi 



The inverse matrix W is then composed of the row vectors Vi, where 

vt = -Niki ( -b n ^ n .fcTZfe-j ^n.fc-3fc-) 

V j^ 1 / 



(54) 



(55) 



One of the eigenvalues will be real and the other two will be a complex conjugate pair: 
k 2 = k\ and k 3 = k%. 

In order to progress further, we need to calculate the eigenvalues. Since we are dealing 
with a 3 x 3 matrix, this is somewhat awkward algebraically. Fortunately, the small j3 
limit is of great interest. The numerical details of the small (3 limit will be given in a later 
subsection. In this limit it is easy to show that the eigenvalues are approximately given by 



k 1 = k* 2 = —D + + 0(f3 2 ), 



(3 2 D 
D 2 + \< 



+ 0((3 4 ). 



(56) 



Since we are still in the regime where the decoherence parameter D is significant, it follows 
that 



Furthermore, 

J>^ dt ' = l + 0((3 2 ). 
Combining these results with Eqs. ( p2[|56] ) we obtain that 



Py(t) 



(3D 



D 2 + A 2 



P z (t) + 0((3 3 



(57) 



(58) 



(59) 



For the antiparticle system we obtain a similar expression, so the neutrino asymmetry evo- 
lution equation (|36|) becomes 



dL a 
dt 



2ru 



2 



D{N a - N a ) D{N a - N s 



D 2 + X 2 



D 2 + X 2 



dp, 



(60) 



14 



having used 



N a -N s 
N Q + N S 



(61) 



and the corresponding antineutrino expression. From Eq.(|60D we see that lepton number is 
constant if the number density distributions of active and sterile species are equal. This is 
to be expected, because oscillations effectively do nothing in that case. 

Since we are working within the R = R = approximation, we now substitute N eq (p, /j,) 
for N a and N eq (p, — /i) for N a . (The chemical potentials of u a and V a are equal and opposite 
because of the rapid Vap a *-> e + e~ process.) Since we are taking the initial lepton number 
to be small, we will also expand N eq (p, fi) to first order in 



to obtain 



u 6n y L a 
T ~ T 3 



N eq {p, n) ~ N eq (p, 0) + N eq (p, 0): 



oP/T 



l + eP/ T T' 

This allows us to form the quantities iV+ and N~ , which are given by 



(62) 



(63) 



N~ 



-(N a + N a ) 



2 ) 



The asymmetry evolution equation 



iV e >,0) + O(L 

12C(3) e vtT 
La n 2 i + eP /T 

) becomes, 



7T- 



s 2 T a(c — b) 



dt 2C(3)T 3 7 [x+ (c-b + a) 2 ][x + (c-b-a) 2 
where c = cos 29 , s = sin2^ ? 



N° q (p,0) + O(Ll). 



■(N+-N+)dp + A + 5, 



(64) 



(65) 



x 





2 


X = 




Am 2 


) 


Am 2 



and 



(r a ) 



p 

(P)o 



(66) 



(67) 



is the collision rate with the chemical potential set to zero. The quantity A is an 0(L a ) 
correction term given by 



A 



7T 



s 2 r [x + a 2 + (b-c) 



4C(3)T 3 J [x + (c - b + a) 2 ] [x + (c-b- a) 2 ] 



(N--N;)d P , 



(68) 



and S is an additional 0(L a ) correction term which arises from allowing D to be different 
from D, 
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4C(3)T 3 J [x + (c - b + a) 2 } [x + (c-b- a) 2 ] 

where xq is x as in Eq. flBlf ) but with r in place of Y a . The evolution equation (J^), including 
the correction term A but excluding 5, was studied in detail in Ref. 0. From Eq.(|65|) we 
see that the effect of having different collision rates for the neutrinos and antineutrinos is to 
introduce the additional correction, 5, as well as a small correction through x and x to the 
denominators. 

The A and 5 correction terms, as well as the 0(L a ) corrections to x and x in the 
denominators may readily be incorporated into the numerical study of the evolution equation 
(|65|), however to a good approximation we may neglect these corrections. It turns out that 
S only amounts to a small contribution such that 

IAI 



0.1. 



(70) 



The corrections to the denominators are very small, even at a resonance. At the initial 



• ' n 



0(L 



M SW resonance where a ~ 0, (b — c) ~ 0, [x + (c — b + a) 2 } \x + (c — b — a) 2 ] 

The approximate evolution equation ( |65|) has proven to be very useful for gaining an- 
alytical insight and indeed for performing numerical work. The 5 = version was studied 
in depth in Ref. J7|. The derivation presented above shows for the first time that the full 
glory of Eq.(^) follows from first principles (that is, the QKEs). In Ref. f7j, a heuristic 
derivation based on a Pauli-Boltzmann-like approximation was used. A minor issue as- 
sociated with the definition of x arises on closer examination. In Ref. 0, this quantity 
contains an additional term sin 2 29o that is absent in the present definition. It turns out 
that keeping higher order terms in f3 in the derivation of Eq . (|60"D does not in fact reproduce 
this missing term as one might hope. The mystery of the missing sin 2 29 seems to require 
an explanation at a deeper level, as yet unprobed. Note also that when the sterile neu- 
trino number density is negligible, Eq . (|65|) becomes a self-contained non-linear differential 
equation describing the evolution of the single momentum- independent variable L a . It is a 
considerable simplification over the full QKEs which are coupled differential equations for 
the eight momentum-dependent functions P, Po, P and P . 

In order to integrate Eq.(^) when the sterile neutrino number density is not negligible, 
an evolution equation for N s is needed. Using 



and the QKEs with the same approximations as above, we easily obtain 



(71) 



d ( N s (p) ' 
dt \N e( i(p,0) 



--pP (p)P y (p) 

'N a ( P )-N s ( P y 



r, 



N e( i(p, 0) J x + (c-b + ay 
T,s 2 



x 



Ax+ (c-b + a) 2 
N s (p) 



1 - 



N e i(p, 0) 



oP/T 



T 3 1 + eP/ T 



- Zn 



1 - 



N s (p) 
N e *(p, 0) 



(72) 
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as the required evolution equation. The sterile antineutrino equation obviously has a similar 
form with the substitutions a — > —a, L a — > —L a and N s — > N s . Equations (|65|), ( 72|) plus 
the sterile antineutrino equation are a useful set of coupled equations that approximately 
describe the evolution of lepton number. 



C. Qualitative features of lepton number evolution 

For completeness, we now briefly recall the qualitative features of lepton number evolu- 
tion which are revealed by Eq.([65f). For further discussion see Ref. 0. Initially, neutrino 
asymmetries are small and hence the term a is small. A critical factor in the subsequent 
evolution of L a is the overall sign of its derivative. (Note also the distinction between L a 
and L^ a \) If it is of the opposite sign to L^ a \ then the asymmetry L a will evolve such that 

L (a) -> 0. (73) 

If, on the other hand, the derivative is of the same sign as L^ a \ then the evolution of L a 
will be such as to increase L^ a \ For Am 2 > 0, the former situation always obtains. No 
explosive creation of lepton number can occur. However, when Am 2 < the situation is 
more complicated. (Observe that for small mixing, Am 2 < means loosely that the sterile 
neutrino is less massive than the active neutrino). In this case the quantity a has the same 
sign as L^ a \ The quantity b is always positive, but it is a decreasing function of time. The 
sign of the derivative is therefore controlled by (c — b). At a given temperature T, the 
function (c — b) equals zero for neutrinos and antineutrinos of momentum p c , where 



7TM W 

Pc 



Am 2 ( 



2T 2 \ V2((3)A a G F 



(74) 



Neutrinos and antineutrinos with momenta less than p c make a positive contribution to 
the righthand side of Eq. (|65"D , whereas those with momenta greater than p c make a negative 
contribution. The critical momentum p c is very small at high T but increases with decreasing 
T, which means that the derivative changes sign from negative to positive once T is small 
enough for p c to be near the peak of the Fermi-Dirac distribution. (For maximal mixing, 
c = and the derivative is always negative. The creation of lepton number occurs due to 
small angle oscillations.) Since a is initially very small, the equation 

c-b = (75) 

is also the MSW resonance condition for both neutrinos and antineutrinos. The denominator 
in the integrand of Eq . (|65|) is thus at a minimum for neutrinos and antineutrinos with 
p = p c ~ p rcs . The resonance momentum p res is generally found by solving 

c-b±a = 0, (76) 

where the plus (minus) sign pertains to neutrinos (antineutrinos). Provided the vacuum 
mixing angle is small enough, the critical temperature T c at which the sign of the derivative 
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changes is well approximated by taking p c ~ 2.2T C which corresponds to the peak of the 
Fermi-Dirac distribution. From Eq. (|74|) , it is given approximately by 



VrM, 



w 



Am 2 c 



4.4 \ y/2((3)A a G F 



1/3 

/ I Am 2 

~15(18)lcos2#o^^] MeV, (77) 




for v e — v s {y^ T — Vs) oscillations. (Numerically, the critical temperature is found to be 
slightly higher than this estimate |7j. Also, when the vacuum mixing angle is sufficiently 
large, a non-negligible sterile neutrino number density is produced, causing the actual T c to 
differ from this estimate.) 

Consider now the parameter space region 

costf ~ 1, |Am 2 | > 1(T 4 eV 2 , (78) 

chosen so that the mixing angle is small and so that T c occurs above the neutrino decoupling 
temperature. At T ~ T c , explosive exponential growth of lepton number begins^ because 
Eq.flBSD is then of the form 

with the proportionality factor augmented by the MSW resonances in the denominator of 
the integrand. After a short time, lepton number has grown sufficiently for the quantity 
a to no longer be negligible. Non-linearity through the denominator then alters the qual- 
itative character of the evolution from exponential growth into a slower phase governed 
approximately by 

a({p)) ~ ± cos 2^, (80) 

depending on whether the sign of the asymmetry is positive or negative. This equation is 
the resonance condition for neutrinos (antineutrinos) when \a\ ^> \b\. 

As the temperature continues to decrease, the collision dominated epoch begins to give 
way to an epoch in which coherent oscillations take over as the dominant driver of neutrino 
evolution. By the time the neutrino decoupling temperature is reached, coherent oscillations 
certainly dominate over collisions. This is the regime relevant for the BBN epoch. Equation 
(|55|) no longer provides a good description of lepton number evolution. The evolution of lep- 
ton number is, of course, in principle obtained by solving Eqs.(^) and (|8]), using the general 
expression given in Eq . (|30|) for the function R. (Numerically, setting D, b ~ is a good ap- 
proximation in this regime.) This in-principle procedure has, to our knowledge, not yet been 
attempted in practice. Instead, the low temperature calculations that have actually been 
performed have utilised the usual adiabatic approximation for fully coherent matter-affected 



7 Note that while the mixing angle should be small, it must be large enough so that sin 2 29 ~ 
5 x 10~ 10 (eV 2 /| Am 2 |) 1//6 is satisfied, otherwise the partially incoherent oscillations are not strong 
enough to create lepton number. See Ref. |?j for more details. 
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oscillations, and the growth of the asymmetry has been calculated using an equation which 
determines how fast the resonance momentum moves through the distribution converting ac- 
tive into sterile neutrinos. The derivation of this equation from the QKEs will be presented, 
for the first time, in Subsection [II E| . Furthermore, the high-temperature expression for R 
as given in Eq. ( |B4|) has been adopted as an approximation to the much more complicated 
and rigorous expression defined in Eq. (|30|) . We may estimate the theoretical error that this 
approximate repopulation method introduces into the final result by comparing the value of 
the asymmetry obtained in the manner described above with that obtained by integrating 
over the neutrino momentum distributions to find the total number densities. If everything 
is consistent, these should of course agree. However, for the calculations presented in Ref. 
[ToU , we find discrepancies of up to 10 to 20%. The repopulation should not change the total 



asymmetry, only distribute it from the particular momentum state at which it is created, 
across the momentum distribution. The discrepancy, which must be due to inconsistencies 
in the repopulation process, is a good indication of the theoretical error in all of the calcu- 
lations HJT0[| performed thus far. This theoretical error is small enough for the conclusions 
so far reached to be essentially unchanged. 



D. Region of applicability for the approximations 



The validity of Eq. (^) rests on the assumptions that the quantity f3 is small relative to 
the size of D and A, and that their rates of change with respect to time are negligible. We 
now examine the regions where these conditions are satisfied. 

The small f3 approximation is sufficiently accurate provided that 







VD 2 + A 2 



< 1. 



(81) 



The "small (3 expansion" is more properly to be thought of as an expansion in the dimension- 
less quantity (3 / \/D 2 + A 2 . In order to explore the bound, we will replace p by its thermal 
average ~ 3.15T. 

Away from resonance, |A| D, so the small (3 expansion holds provided that \f3\ <C |A|. 
For T > T c , A ~ Am 2 b/2p, which means that the small (3 limit requires 



| Am 2 
eV 2 



T \ ^ 



■sin 2^o < 2 x 10- 7 (6 x 10~ 8 ) (— ) 



(82) 



for a = e and /i/r, respectively. We would like the small (3 limit to hold just prior to lepton 
number creation at T = T c . Putting T = T c , as given by Eq. (|77|) , in Eq.(|82]) we see that 
| Am 2 1 cancels out of the inequality, leaving 



tan 29a -C 1 



(83) 



as the constraint. Note that the use of Eq.([T7|) presupposes that the mixing angle is suf- 
ficiently small for the sterile neutrino number density to be small, so this is not a serious 
constraint. 
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In the centre of the lepton-number creating resonance at T = T c , A momentarily goes to 
zero. The small j3 expansion will be valid in the centre of the resonance, provided the more 
stringent constraint |/?| <C D holds. This translates into the requirement that 



| Am 2 
eV 2 



sin 26 <4 x 10" iU y, 



MeV 



10 



-9 



T c \ 6 
MeV/ 



(84) 



Using Eq.([77]), we find that |Am 2 | cancels out of the inequality, leaving the constraint 

tan2# <10~ 2 (85) 



on the vacuum mixing angle. 

We now discuss the applicability of the adiabatic-like approximation which assumes the 
complete negligibility of the term that appears in Eq. fl48|) . Evaluated explicitly in 

the small (3 limit, this hitherto discarded quantity takes on the form 



U 



.dlC 

at 





( w 


Y* 


-Z* 




Y 


W* 


-z 




\ Z* 


z 


X 



(86) 



with 
W 



1 (3D D — i3A dp 



+ 



X 



Y 



2XD 2 

i(3 2 J 
D 2 + A 2 I 

-2(3D 2 



X 2 X + iD dt 
{ dD 



dt 

d(3 



J_ (X-iD)(D 2 - 
4A + (D 2 + A 2 ) 2 



A 2 ) 



+ 



dX 
dt 



D 2XD(X 
+ 



iD) 



4A 



+ 



2(3 2 D 



(D 2 + A 2 ) 2 dt (D 2 + A 2 ) 3 
i (3D d{3 i (3 2 



dD 



dX 



2X D 2 + A 2 dt 
1 1 d(3 
71 (A -iD) ~dl ~ ^/2{X-iD) 



(D - X )—— + 2XD— 

dt 

+ 0((3 3 



dt 
dX^ 



(D 2 
0(P 3 ) 



+ A 2 ) 2 



f dD 

4A 2 D 2 + A 2 [ ~dt + dt 



P 



.dD dX\ ^, n2 



0(P 3 ), 



(87) 



to the lowest order in p. These non- vanishing diagonal and off-diagonal entries contribute 
to corrections to both the eigenvalues and eigenvectors of the matrix /Q in Eq. (T47|) . Thus 
our immediate task is to ensure that these corrections are small in the region of interest. 

The 33 entry of the matrix UdlA~ x jdt represents a "first order" correction to the eigen- 
value kz- Imposing the condition 



h 

and rendering it into a more illuminating form 

1 dp 



< 1, 



(88) 



(D + A ) » 



2D p^ {d2 + x2) - 2 



, 9 , 9n dD , dX 
(D - X )—— + 2XD — 
v 1 dt dt 



(89) 
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we find that 



T > 



4.3 



1/3 



MeV ~ 3 MeV, 



(90) 



where the momentum p has been replaced with the mean momentum (p). Note that in 
deriving the above, we have set A ~ 0, i.e., the resonance condition, and we have used the 
relation t ~ mpi anc k/llT 2 . Thus Eq.(pO|) purports the failure of the static approximation for 
resonances occurring at temperatures below roughly 3 MeV. Together with the resonance 
condition c ~ b ~ 1 , Eq. (|90|) effectively pi lower limit on |Am 2 |, that is 



| Am 2 
eV 2 



> 8.7 x 1Q- 6 A Q ~ 10~ 4 (4 x 10" 



(91) 



for v P 



- v s 



oscillations. 



Additional bounds arise from evaluating Eq.(p9|) alternatively at A ~ D, i.e., by requiring 



that 



D 2 » 



D_d£ 
~]3~dt 



~dt 



(92) 



The relative importance of the two terms above is not immediately obvious by inspection. 
Suppose the first term dominates in the region of interest. The resulting inequality is but 
a less severe version of Eq. (p0|) . Supposing now that the second term is predominant and 
that there are no accidental cancellations such that the inequality roughly reduces to two 
separate conditions 



D > 



D > 



11T 2 



fTlpianck 
11T 3 



D 

2 + 
Am 2 



Am 2 



2p 
a dL^ 



(3b ± 2a) 



2m Planck 2p LW dT 



(93) 



While the former is essentially another encrypted bound on the temperatures at which our 
approximations are valid, similar to that obtained earlier in Eq.(^), the latter amounts to 
demanding that 



dT 



« 3 x 10- 12 y 2 



T 



MeV J MeV 



(94) 



at the mean momentum. This condition agrees remarkably well with that obtained in Ref. 
|0] and should be checked when integrating the static approximation equation (p5|) for self- 
consistency. 

We shall not reproduce here the constraints resulting from the (UdU^ 1 / dt)u and 
{UdU^ 1 / dt) 22 corrections to the eigenvalues k\ and k 2 respectively, as they appear to be 
much less stringent than those arising from earlier considerations. Indeed, this is to be 
expected given that the correction terms are only of the order f3 2 where the small (3 limit is 
effective. 
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Turning our attention now to the effects of U.dlA~ x jdt on the instantaneous eigenvectors, 
we observe that, using a perturbative analysis, the "first order" corrections to Q[ (t), Q 2 °\t) 
and Qs(t) take on the forms 



Q { i\t) 
Q?(t) 



U 



k\ - k 2 



U 



au- 1 
dt 



12^,(0) 



k 2 - h 



Qnt) + 



k\ - k 3 



Q 3 (t), 



("ST 1 ), 



k 2 - h 

-l 



k 3 - ki 



QT } (t) + 



23^(0) 



fa - ko 



V2 



Q2 (t) 



exp(#)Q! 0) (t) + exp(-i0)Q^(t) 



(0), 



(95) 



where 
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exp(2<; 



fei 



12. 



(96) 



and we have neglected terms proportional to (3 2 . Thus the single requirement that must be 
satisfied in this case in the small (3 limit is 
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<c 1. 



(97) 



It turns out that the above condition does not in fact lead to more severe bounds than those 
following from other requirements. 



E. Oscillation dominated asymmetry evolution 

At low temperatures, such that D <C \/3\, |A|, the collision rate becomes much slower and 
the generation of lepton number is dominated by oscillations. Taking the D = limit, we 
can demonstrate that the QKEs reduce to usual expression for MSW transitions. 

We begin by setting D = in Eq. (|44j) . In this case, setting dU^ 1 /dt = exactly 



corresponds to the usual adiabatic approximation for MSW evolution. The eigenvalues now 
become 



k l = k* 2 = i^f3 2 + A 2 , k 3 = 0, (98) 
and the matrix U^ 1 {= U^) is given by 
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u 



-1 



A A y/20' 

-i^(3 2 + X 2 i^P 1 + A 2 
-(3 -(3 V2X 



'2(/3 2 + A 2 ) 

With the initial conditions P x (0) ~ 0, P y (0) — 0, the solution for P z (t) is given by 

P z (t) ~ ^W^(t) exp [ T fe^)^] W l7 (0)P 7 (0) 



(99) 



1,7 



A(f) 



A(0) 



where C26» ro = cos20 m and 



C2eJt)c 2drn (0)P z (0), 



sin 2 (2^ 



s 2 + (c-6 + a) 2 ' 
The number density of sterile neutrinos is then given by 



1 



N 8 (t) ^ N a (0)~ [1 - c 2dm (t)c 2em (0)] + N s (0) 1 - - [1 - c 2 e m (t)c 2em (0)] . 



1 



(100) 



(101) 



(102) 



Since the adiabatic transition probability for a neutrino propagating through a medium 
where the mixing angle changes from 9 m (0) to 9 m (t) is given by 



Prob(f 



1 



[1 - c 2 e rn {t)c 2 e n 



we have 



N s (t) ~ iV a (0)Prob(z/ a -> v s ) + iV s (0)Prob(z/ s 



(103) 



(104) 



so, as expected, the adiabatic approximation for the QKEs is equivalent to the adiabatic 
limit of the MSW effect if we turn off the decohering collisions. 

The low temperature equation, used in previous work PJlOl to calculate the growth 
of asymmetry by determining how quickly the resonance momentum moves through the 
distribution, can be derived from the QKEs in the D = limit. 

Since we have b ~ 0, the resonance condition c ± a ~ can only be attained for either 
the neutrinos or antineutrinos, depending on the sign of the asymmetry. Assuming for 
defmiteness that L a > 0, which implies a > 0, only the antineutrino resonance will contribute 
to the growth of lepton number, so that 



dT dT n 



4^— i d p 



~P Q (1 + P Z )N C * - ip (l + P z )N e * 



dpN^P dP z 
2n 7 dT 



(105) 



where ~ has been assumed. From Eq.( |100| ) we find 
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iF ' - PM L.f.-'fl.J f o . r. * n ) ^ (¥ ) doe) 



dT zy ' \[s 2 + (c - a(0)) 2 }V 2 J \[s 2 + (c - a(t)) 2 } 3 / 2 J p res dT \ T 
where we have used a/c = p/p res and = —j^-j^ (^-p)- This results in 
dL a T d fp 



dT 2n 7 p res dT \ T 

x CZ" dpN °™ {vTjf^kw) { p + A)W> ) • (107) 

where 5p res is the resonance width, since the last term in Eq. (|107|) ensures the integral is 
strongly peaked about the resonance momentum. Averaging the integrand over the width 
of the resonance leads to 

dL a _ N a (p ms + Sp rcs ) - N s (p ves + 6p ieB ) T d /p res \ 



dT ra 7 dT \ T J s 2 + (c5p ICS /p vcs ) 2 

so that in the limit of zero resonance width we recover 

where the N a — N s term is interpreted as evaluated just before the resonance passes a 
particular momentum value. This term expresses the difference between the number of 
V a and V s which are converted at the resonance, while the factor ^ {^^j accounts for how 
quickly the resonance momentum moves through the distribution converting active to sterile 
neutrinos. See Refs. P,|TU[ for applications of this type of equation. The derivation given 
above for Eq.( |109| ) is the first justification of it from first principles. 



III. MIRROR NEUTRINOS INSTEAD OF STERILE NEUTRINOS 

We now briefly discuss how the foregoing must be altered to deal with mirror instead 
of strictly sterile neutrinos. Mirror neutrinos form the neutral lepton sector of the mirror 
fields postulated in the Exact Parity Model (EPM) ||20|| . The EPM restores exact parity 
invariance to the microscopic world by parity-doubling the Standard Model. Interestingly, 
if neutrinos and mirror neutrinos have mass, and if the two sector mix, then in the absence 
of intergenerational mixing the mass eigenstates are maximal superpositions of ordinary 
and mirror neutrinos. Maximal oscillations of with its mirror partner 1/ can explain 
the atmospheric neutrino anomaly, while maximal oscillation of u e with its mirror partner 
u' e can similarly solve the solar neutrino problem. The maximal mixing of v a with v' a is a 
consequence of the unbroken parity symmetry. This form of parity symmetry is arguably 
one of the most credible explanations for the maximal mixing of muon-neutrinos observed in 
atmospheric neutrino experiments. The LSND anomaly can be accomodated by switching 
on small intergenerational mixing with the appropriate mass hierarchy. 

The issue of whether or not the EPM explanation of the neutrino anomalies is consistent 
with Big Bang Nucleosynthesis is an important one. We assume that in the very early 
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stages of the Big Bang ordinary matter predominates over mirror matter. (See Refs. |2T| for 
speculations about how this could arise.) The question is then whether interactions between 
the ordinary and mirror sectors at later times overcreates mirror matter, spoiling BBN. 
Focussing on ordinary-mirror neutrino mixing, we note that a sufficiently large (~ 1CT 5 ) pre- 
existing neutrino asymmetry can always be invoked to suppress ordinary-mirror oscillations, 
thus leading to standard BBN despite the presence of the mirror sector ||. However, it is 
more interesting to suppose that the neutrino sector of the EPM "saves itself" through the 
production of large neutrino asymmetries in the manner discussed above. In order to do 
this analysis, the QKEs have to be modified to take into account the interactions of mirror 
neutrinos amongst themselves and with other mirror species such as the mirror photon and 
the mirror electron. See Ref. || for the proof that the EPM is consistent with BBN in an 
interesting region of parameter space. 

Mirror neutrinos feel mirror weak interactions. The interaction strength is given by Gp, 
since exact parity symmetry imposes equality of coupling constants and gauge boson masses 
between the sectors. We may classify the additional terms in the QKEs through powers of 
Gp: coherent forward scattering induces 0(Gp) terms, while incoherent scattering induces 
terms of 0(G 2 F ) and higher. Provided that the number densities of mirror neutrinos remains 
low, all terms other than the 0(Gp) coherent forward scattering effects can be neglected. If 
the number densities of mirror neutrinos become appreciable, then many of the complexities 
inherent in active-active scattering become an issue. 

So, provided one is working in a parameter space region where the mirror neutrino 
number densities remain small, it is consistent to alter the QKEs of the strictly sterile 
neutrino case by the simple substitution, 

£(«)_> £ (a >-Z/<fl, (110) 

for v a — v'p oscillations (a,/3 = e,/i, r). The quantity L'^' is simply the mirror version of 
L^>. This substitution fully incorporates the coherent forward scattering of mirror neutrinos 
off mirror neutrinos induced by mirror weak interactions. 

It is important to understand why these 0(Gp) mirror neutrino self- interactions are as 
important as 0(Gp) neutrino self-interactions. The term a(p) in the effective matter poten- 
tial is equal to the sum of terms directly proportional to the difference in number densities 
between neutrinos and antineutrinos, and the difference in number densities between mirror 
neutrinos and mirror antineutrinos. These number density differences will be of the same 
order of magnitude simply because mirror neutrinos and antineutrinos are produced, via 
partially incoherent oscillations, from ordinary neutrinos and antineutrinos. Since we know 
that the a(p) term dominates the b(p) term once significant neutrino asymmetries have been 
created, and since coherent effects anyway become more and more important as the temper- 
ature approaches the BBN epoch, the 0(Gf) mirror neutrino self- interactions are of critical 
importance. The fact that they are completely non-negligible distinguishes the cosmology 
of mirror neutrinos from that of strictly sterile neutrinos p2|. 
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IV. MULTI-FLAVOUR OSCILLATIONS 



The purpose of this section is to begin the study of multi-flavour effects. In particular, we 
wish to determine the conditions under which it is appropriate to break up a multi-flavour 
situation into effective two-flavour subsystems, where two-flavour expressions can be applied. 
An understanding of when multi-flavour effects may arise is important, because while two- 
flavour subsystems are a useful approximation, a realistic situation will always involve three 
or more flavours. A multi-flavour system may involve features that do not occur in simpler 
two-flavour systems. We would expect, for instance, genuine multi-flavour effects to arise in 
systems in which there are "overlapping" resonances. With a general multi-flavour system, 
there is also the possibility of explicit CP violation due to complex phases in the neutrino 
mixing matrix. 

Partially incoherent oscillations in a multi-flavour system have not previously been stud- 
ied from first principles. In all prior work on lepton asymmetries, multi-flavour systems were 
dealt with in terms of two-flavour subsystems. 

The system we shall consider here, chosen both for definiteness and importance, consists 
of maximally mixed a sterile neutrino u s , and a heavier v T which has small mixing with 
Up and u s . This particular system is well motivated by the atmospheric neutrino anomaly, 
with the parameter space region sin2# MS ~ I and Am^ s ~ 10 -2 — KR 3 eV 2 leading to a 
resolution of the anomaly via — > v s oscillations. The v T ,^,s system was studied in Ref. 
|7]], where the role of the heavier tau neutrino, taken to have a mass in the range where 
Am 2 s = f — fOOOeV 2 , is to generate an L T asymmetry via oscillations with the sterile neutrino 
as per the dynamics discussed in the previous Section. The L T asymmetry generates a large 
Wolfenstein term in the effective potential for the — v s subsystem, which then acts to 
suppress oscillations between the muon and sterile neutrinos. For a range of parameters, this 
prevents the maximal — v s oscillations from bringing the sterile neutrinos into equilibrium, 
thus circumventing the supposedly stringent BBN bounds on the mixing of sterile and active 
neutrinos. This phenomenon is vital for reconciling the z/ M — > v s solution to the atmospheric 
neutrino problem with cosmology. Note, however, that in previous analyses the u T — v s and 

— v s resonances were dealt with in a two-flavour manner; possible three-flavour effects 
were ignored. Given that the resonances are well separated, this is expected to be a good 
approximation. However, in order to verify this is the case, we need to study it within a full 
three-flavour framework. 

We begin by writing down the evolution equations for a three-flavour partially incoher- 
ent system using the density matrix formalism. While it is in a sense straightforward to 
generalise the Quantum Kinetic Equations to a three-flavour system, qualitatively different 
complexities are encountered due to the presence of active-active neutrino oscillations in 
addition to the more readily understood active-sterile oscillations. These active-active com- 
plications result in equations from which it is difficult to obtain analytical insight .0 However, 
under the momentum averaged approximation the equations simplify significantly, allowing 



8 See Ref. [ 14 1 for details of the (momentum dependent) Quantum Kinetic Equations for active- 



active oscillations. 
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progress to be made. Although going to the momentum averaged approximation will in- 
troduce some errors, it is probably a reasonable approach. In any case, we adopt it here 
because it makes possible a first attempt at a first-principles analysis of a three-flavour 
system. (The differences between including and excluding the momentum spread for two- 
flavour active-sterile oscillations were discussed at length in Ref. J/J. The qualitative features 
of the evolution were found to be preserved in the mean momentum approximation.) 

The three-flavour momentum averaged density matrix (p) is parameterised in terms of 



the SU(3) Gell-Mann matrices Aj 



8 such that 



(p) 
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The neutrino number densities are related to the diagonal entries of (p): 
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The evolution of a momentum averaged density matrix is described by the Quantum Rate 
Equations (QREs) rather than the Quantum Kinetic Equations. Generalising the two- 



flavour QREs derived in Ref. [14], which involves the lengthy but straightforward procedure 
of substituting Eq. ( |1 1 2j ) into the general QKE expression given by Eq.(22) of Ref. |H, we 
find the three-flavour QREs describing the considered system to be, 



dt 



V x P - D(P lX ! + P 2 x 2 + P e x 6 + P 7 x 7 ) - D'(P 4 x 4 + P 5 x 5 ) 
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;ii4) 



with the evolution of Pq, which is related to the total number density of z/ r , and u s , given 
by 



dPo 2 



(115) 
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where R a is the repopulation function |L4 
1 



J2 ( r (^' U a"a)) 



hjUjUj - n Va n Va 



2 S G * 



P4P4 - ft? 



B r 5 



;ii6) 



where h e = \, h v — 1 and the values for the collision rates (T(jj — > VqVo)) can be found 
in Refs. [ff3,B31. The n's are normalised such that nf q = n eq and n*; q = 2n eq . The G terms 
arise from the mixed active-active part of the density matrix, with G e = 0.26Gf/T 5 and 
G Ve = 0.51G 2 F T 5 . 

The quantity D = |(r a ) is the damping parameter, with T a as given by Eq. fl27|) , while 
D' is the equivalent parameter for the active-active v T — oscillations: 



D' ~ 1.2GlT 5 . 



(117) 



The parameter C in Eq.( |114[ ) is a damping-like term which couples the neutrino and an- 
tineutrino density matrices through the off-diagonal p^ T and p-^ elements, and only arises 



in the case of mixed active-active neutrinos. For v T — oscillations W\ 



C ~ 1.8GIT 5 . 



(118) 



The quantity H is given by 



H 



n eq 



dk'dp'dpdk5E(k + p — k' — p')f(k)f(p) 



xE[^W*).^(p)ljX^)J(^))][^OX*0J(i/)k(*)^r(p))], 



(119) 



where / dp = / c? 3 p. This is an order G F T 5 quantity, however it will be small compared 
to D and C, provided that is small. 
We define 



(V x P) fc = T^pi/y*, 
where / lJ are the SU(3) structure constants. Finally, V is given by 

V = 2ReE TS x 1 - 2lmE TS x 2 + (E TT - E ss )x 3 + 2ReE T ^x A - 2lmE Tfl x, 

+ 2ReE s »x 6 - 2lmE s ^x 7 + ^{E TT + E ss - 2E^)x 8 , 

v3 

with 

where u al3 are the (vacuum) energy eigenvalues in the flavour basis, 



(120) 



CO 



a/3 



2p 



f/diag(m^, m 2 ,, ml)W , 



;i2i) 



(122) 



(123) 
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with U being the mixing matrix relating the neutrino mass and flavour eigenstates. V alS are 
the effective potential terms which, as well as the usual diagonal potential terms (V aa = 
a + b)), also include off-diagonal contributions given by 



Am 
2p 



yttr 



V2G 



(P M r - Pw) 



14C(4) 2T 



3M 2 ' Ptir 



(124) 



The appearance of off-diagonal terms in the effective potential is a consequence of active- 
active oscillations. They depend on the quantities P4 and P5 which parameterise the mixed 
v T and i/p states in the density matrix. Note, however, that the number densities of v T and 
Up are equal (to order L). As we shall see, this results in P 4 and P 5 being suppressed, with 
respect to Pi & P2 and Pq & P7 which parameterise the mixed u Tj8 and v^ s subsystems, 
respectively. 

For simplicity, the v T — mixing angle will be set to zero. Observe, however, that 
there will still be a small effective mixing between v T and indirectly through v s . We 
parameterise the mixing matrix U as 



U 




(125) 



where <fi is the u T — v s mixing angle, and the — v s mixing angle has been fixed at ?. 

Unlike the two-flavour case, the neutrino and anti-neutrino density matrices are coupled 
through the C term in Eq.( |114j ), so that we must consider a system of 16 coupled differential 
equations for P 1; . . . , P 8 , P 1; . . . , P 8 . As in the two-flavour case, we adopt the ^ ~ 
approximation. This allows Eq.( |114|) to be expressed in the form: 
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/CP, 
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where P 



(Pi, ...P 8 , Pi, 



Pg) and /C is the matrix given by 
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where M. and C are the submatrices 



M 



( 


-D 


-A 








b 


\5i — Hr 


\5r — Hi 





\ 




A 


—D 


-P 


-b 





\5r — Hi 


-bi+ H R 

















bi 


bn 





-b 












b 


-\5i 


-D' 


—0 





-b 








-b 





-bn 


a 


—D' 


b 













— — Hi 








-b 


-D 


— e 









—\$r — Hi 


\5l + Hr 


b 


b 





e 


-D 


2 ' 




\ 













—Or 





2 ' 





/ 



;i28) 



29 



c 



/O 





c 







\0 



\ 







-C 





/ 



(129) 



and A4 is obtained from M. via the replacement A — > A, o 
and H —> H. We use the notation, 5j = Im<5, 5 R = ReS, Hj 

Am 2 



W, e -f e, D -> D, £>' - 
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where Am 2 = Am 2 and r 



Am 2 J Am 2 . s 



< 1. 



, s/ ... .... The parameters /3 and A take roughly 

the same form as for the two-flavour active-sterile system, and are related to the v T — v s 
oscillations. The quantities 5 and 7 are the equivalent of j3 for the v T — and z/ M — v s 
oscillations, respectively, while a and e are the v T — and — v s analogues of A. 

From Eq.( |114| ) and the approximations adopted as per the discussion above, the evolution 
of the lepton numbers L T and are found to be 



cq 



dt 2n„ 



P (PP 2 + 6 I P A + 5 R P 5 ) - P ((3P 2 + +5jP, + 5 R P 5 



~dT 



2n„ 



-P ^P 7 + +61P4 + S R P 5 ) + P (^P 7 + +S1P4 + 5 R P 5 



;i3i) 



To solve these equations, P 2 , -P4, P5 and P-j must be determined as functions of time. 

A first inspection of Eqs. (|127|) , ( |128|) and (|129|) gives the impression that the various 



oscillation modes are coupled together in a non-trivial manner. However, under the assump- 
tions made, the solution reveals that the three oscillation modes effectively decouple to first 
order in the small parameters (3, 5, 7 and H. 

To solve these sixteen coupled equations, the matrix will be diagonalised in the small 
/3,5,"/,H limit, in analogy to the two-flavour case. This is done by treating the /3, 5, 7, if 
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terms as a perturbation to the zeroth order solution obtained with {3 = 5 = ^ = 11 = 0. The 
parameters (3 and 7 may be of approximately the same order of magnitude, depending on 
the mass-squared differences and mixing angles. The size of 5 depends on the size of the off- 
diagonal potential term V T ^ which, as with H, is dependent on the values of P4 and P5. Our 
perturbative solution relies on the assumption that P4 and P5 are small. This is expected 
for the physical reason of decoherent damping through D', except at low temperatures (say 
in the neutrino decoupling regime). The solution obtained with this assumption may then 
be checked for self-consistency. 

For the purposes of the perturbation, we treat (3, 5, 7, H as being of the same order. 
Again, as with the two-flavour case, we adopt the adiabatic-like approximation that U^- ~ 
where IA is the time dependent matrix which instantaneously diagonalises the matrix in 
Eq.( |127|) . Solving for the appropriate eigenvectors to first order in (3, 5, 7, H yields 



{3D 



D 2 + A 2 V P n cq 



P 5 (t) = Y 5 
Pr(t) 



P n ec i 
n,, - n h 



'n Vr - n v 



P n ec i 
7P /n,,.. - //, 



P n ef i 



D 2 + e 2 V P n e ^ 



(132) 



where I4 5, Z 4 5 are order 5 terms and are functions of the parameters D', a, a and C, arising 
from the diagonalisation of the P4, P5, P4, P5 submatrix 
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We see that P5 (and similarly P 4 ) is suppressed not only because 5 is small but also 
since n Vr — n u = O(L). This ensures that V T[1 and H remain small, consistent with the 
assumptions made. The proportionality of P4 and P 5 to n Vlx — ti Ut is one manifestation of 
the unimportance of oscillations between species of roughly equal number densities. 

These expressions in Eq. (|132|) resemble those that would be obtained by breaking the 
system down into two two-flavour subsystems v T + v s and + v s . The evolution equations 
for L T and L M are coupled only via the effective potentials and n Vs . In the initial stages 
of lepton number creation, when the sterile neutrino number density can be neglected, the 
equations are in essence coupled only through the Independent effective potentials. The 
equations for lepton number become 
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with the P4 and P5 terms now neglected. These types of equations, generalised in the 
obvious way to incorporate the thermal momentum distribution, have been used in practical 
calculations |7],§]. The above analysis is the first to fully justify their use on the basis of first 
principles (admittedly in the Quantum Rate Equation approximation). 

The various oscillation modes would seem only to be coupled via higher order terms in 
(3, 7, 5, H. In fact the restrictions we imposed on the mixing angles can be relaxed, and 
the only difference between Eq.( |134j ) and the two-flavour equations used in Refs. MM are 
the small corrections (in terms of the mass-squared differences and mixing angles) to the 
parameters /3, A, 7 and e. 

This two-flavour subsystem solution will not be valid when the resonances of different 
subsystems are close enough to overlap. For example, the v T — v s and z/ M — v s resonances 
begin to overlap if A ~ e. In this case some of the eigenvalues of the matrix in Eq.( |128| ) 
become degenerate and hence the perturbative method used to obtain the solution may 
break down. 

The study of the v T ^ s system, with z/ M — + v s parameters set by the atmospheric neutrino 
anomaly, has recently been re-examined by Foot |TT|] . In this work, the "pairwise two-flavour 
approximation" is used, but with each two-flavour active-sterile subsystem treated using the 
full two-flavour QKEs. 

V. CONCLUSION 

Active-sterile neutrino oscillations will play an important role in early universe cosmol- 
ogy if light sterile neutrinos exist. The combined solar neutrino, atmospheric neutrino and 
LSND anomalies require at least one light sterile flavour. The partially incoherent neu- 
trino oscillations that occur in the early universe are described using the Quantum Kinetic 
Equations. These equations provide an in-principle means of calculating the cosmological 
consequences of light sterile neutrinos, with Big Bang Nucleosynthesis being an important 
concern. However, their qualitative physical consequences are not always easy to extract, 
and they are computationally demanding. Approximation schemes are therefore welcome. 

The main insights of this paper are twofold: 

1. An adiabatic approximation for the partially incoherent oscillations described by the 
QKEs has been developed. This clarifies the origin of the "static approximation" 
discussed in Ref. 0, revealing it to be conceptually related to the usual adiabatic 
approximation of matter-affected oscillations. In the absence of collisions, we have 
explicitly shown that it reduces to the usual adiabatic approximation. When collisions 
are important, we have used it to rederive some very useful approximate evolution 
equations for lepton number in the small (3 limit. The evolution of lepton number is 
of particular importance because of the central role played by the Wolfenstein term 
in the effective matter potential in suppressing sterile neutrino production prior to 
BBN. The systematic approach to approximating the QKEs introduced in this paper 
suggests further development which could improve, in a hopefully practical way, on 
the adiabatic approximation for partially incoherent oscillations. 
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2. A first principles treatment of partially incoherent three- flavour oscillations has been 
attempted for the first time. Focussing on a v T ,^,s subsystem with maximally mixed 
and u s , important for resolving the atmospheric neutrino anomaly, the three-flavour 
Quantum Rate Equations were explicitly written down. The QREs are obtained from 
the QKEs by approximating the evolution of the neutrino momentum distribution 
by the evolution of neutrinos having the mean momentum. By employing a similar 
adiabatic-like approximation in the small /3, 7, 5, H limit (three-flavour static approxi- 
mation) to this case, it was demonstrated that the three-flavour system separated into 
two two-flavour subsystems, v T + v s and + u s , coupled only through the dependence 
on the Wolfenstein term on family lepton numbers, and through the common u s en- 
semble. The "coupled two-flavour subsystem" approach had previously been used in 
studies which concluded that the maximal — > v s solution to the atmospheric neu- 
trino anomaly was consistent with BBN provided that the v T — v s oscillation parameters 
lay in a particular region. The conclusion that there is no cosmological objection to 
the — > v s solution to the atmospheric anomaly can thus be made with even more 
confidence. 

An exciting new era in fundamental physics has begun, ushered in by the beautiful 
atmospheric neutrino results of SuperKamiokande. We await with great interest news from 
the Sudbury Neutrino Observatory |24|, and further data from SuperKamiokande, regarding 



the existence or otherwise of light sterile neutrinos. Should they exist, then a new synergy 
between the microscopic and macroscopic worlds will be revealed through the evolution of 
partially coherent active-sterile (or ordinary- mirror) neutrino oscillations. 
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APPENDIX A: CALCULATION OF DECOHERENCE FUNTION 



In this section the calculation of the decoherence function D(k) is outlined. The deco- 
herence funtion was derived in Ref. [13] as 

D(k) = nJ dk'dp'dpdEik + p-k'-p^^iV'M^^ip^k^^ip'))/^) 

j 

+ V\v a (k),V a (p)\j(k')J(p'))f Va (p)}, (Al) 
where / dp = / d 3 p, and the /'s are momentum distribution functions with 

ffW)= , * v (A2) 
1 + exp(^) 

and f£ q = f% q (p,~p). Note that this expression neglects Pauli blocking factors of (1 — 
/(p'))(l — f(k')). The sum over j in Eq.([Al|) includes all weakly interacting particle species 
in the background plasma. 

The matrix elements V(i(k), j(p)\m(k'),n(p')) are related to weak interaction matrix 
elements M(i(k), j(p)\m(k'),n(p')) via 

V 2 (z(k),j(p)\m(k'),n(p')) = 2k ^ l2p A k +p-k'- p')M 2 (i(k)j(p)\m(k')n(p')). (A3) 

The matrix elements M(ij\mn) may readily be evaluated as they are simply four- Fermi 
interactions, and the integration over k' and p' performed, to obtain 

where the mass of the electron has been neglected, and the Aj's are coefficients which are 
given in terms of sin 2 9w, and can be found, for example, in Refs. |14^3| . 

If we assume that all weakly interacting species are in thermal equilibrium with zero 
chemical potential, then we obtain 

„(*) = (A5, 

which, if we set the momentum k equal to its thermal average (k) Q , reduces to the standard 
expression in terms of the thermally averaged collision rate (D) — |(r a ). 

To allow for a nonzero lepton asymmetry, we will assume that u a and V a may have 
non-zero chemical potential, while all other particle species will be assumed to have zero 
chemical potentials. This results in small correction terms to Eq. ([A5[) 



D(k) = \g 2 f T*^ L + uj^ + v a ^ + w a (f) 2 + x a (^) 2 + O 



T 



(A6) 



where y e ~ 4.0, y^ T ~ 2.9, u a ~ 0.72, v e ~ 1.0, v^ T ~ 0.8, w a ~ 0.33, x e ~ 0.46, and 



x^ ~ 0.36. 
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APPENDIX B: REPOPULATION 



We outline the approximations used in determining the form of the repopulation function 
R(k). The general form of this function was derived in Ref. and is given by 



R{k) 



2tt 



J dk'dp'dpSE(k + p — k' — p') x 



E^ 2 K(A;),j(^)l^(^),ib , ))(^(^)/ib / ) - UMfM) 
j 

+ v 2 (v a (k),v a (p)\j(k , ),j(p , ))(fj(k , )fj(p') - UMhM)]- 



(Bi) 



Note that both inelastic (annihilation) and elastic (scattering) processes are included in the 
expression for R(k), because both contribute to the rate at which a certain momentum state 
is refilled. However, in the momentum averaged limit, the elastic contribution will vanish 
since in that case we are only interested in the total number of particles, and not how they 
are spread across the momentum distribution.^ 

The equation for R(k) is a Pauli-Boltzmann equation, and momentum states are refilled 
in such a way as to drive them toward equilibrium. The term (f(k')f(p') — f{k)f(p)) 7^ if 
all the /'s are given by equilibrium Fermi-Dirac distributions, but this is just a consequence 
of neglecting the Pauli blocking factors. The general form of R(k) cannot be calculated an- 
alytically, and numerically would require significant computing power, so it is useful to look 
at approximations under which it may be simplified. If we assume that all the distributions 
are in thermal equilibrium except for the state which is being refilled [i.e. fv a {k)}, then we 
may replace fj{k')fj(p') in Eq. ([BT]) with f^{k)f^(p). This makes sense in terms of pushing 
everything toward equilibrium, or alternatively, 
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since Pauli blocking factors are neglected. With this approximation, 
2tt 
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so that we obtain, 
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9 The expression given in [14| appears to neglect elastic processes. 
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The general expression for the repopulation function given by Eq. ( [Bl|) is numerically inten- 
sive, so the form given by Eq.(B4) allows a significant simplification to numerical calcula- 
tions. 
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